Efficient entanglement length measurements for photonic cluster state sources 
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We present a procedure for confirming the generation of long strings of photons in an entangled 
(cluster) state that does not rely on complete state tomography and that works even at low collection 
efficiencies. The scheme has the added advantage of being passive - it does not require switching of 
optical elements to perform measurements in different bases for instance. 
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Recently, a scheme for preparing on-demand photonic 
cluster state strings (a "cluster state machine gun" Vwas 
proposed [l| and adapted to various architectures 
Such a device would ideally allow for the production of 
hundreds of entangled photons entangled in a ID cluster 
state, a technology that would dramatically change the 
face of optical quantum information processing. How- 
ever, it immediately brings to the fore the challenge of 
benchmarking and characterizing the output light: How 
docs one verify that the emitted photons are actually en- 
tangled as expected, and how many photons are entan- 
gled before inevitable noise destroys the quantum corre- 
lations? 

The simplest solution would be to perform a full to- 
mography on K photons to deduce the reduced density 
matrix pK of these K photons. From px one can try 
and calculate the entanglement between the photons us- 
ing known entanglement measures. The problem is that 
given current experimental setups, K must be quite small 
because of collection & photodetection efficiencies which 
reduce the number of sequential photons that are mea- 
sured. 

For example, denoting by pd the probability of photon 
collection and detection, the probability of detecting K 
photons in a row is {pd)^ ■ Furthermore, there are 2^^ 
variables in the density matrix of K photons. In order to 
do the full tomography, even given a good passive exper- 
imental setup, one will need at least measurements 
of K photons in a row to get a measurement of every one 
of the variables (actually a lot more, but this will be a 
lower bound). Therefore, we get that the amount of data 
needed for one measurement of each variable is bounded 
below by 
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with N being the number of photons emitted. For one 
photon emitted per nanosecond with collection efficien- 
cies of 10%, 50% and 90% and a measurement time of 10 
seconds, this yields 6, 11, 15 as the best case number of 
entangled photons in the reconstructed cluster state. 

In addition to the issues of collection efficiencies, short 
(nanosecond) time scales between photons are needed for 
the quantum dot example considered in Ref. [l| to achieve 
minimal decoherence. As such allowing active optical 
elements that can be switched between photons in the 
setup will be very difficult. 
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FIG. 1. The proposed experimental setup for measuring the 
entanglement length in a passive manner. The reflectivity of 
the beam-splitters can be adjusted so the probability of mea- 
suring the photon in each of the Pauli bases is controllable. 
The default detector polarization basis is assumed to be the 
horizontal and vertical polarizations. 



In this paper we propose an experimental setup that al- 
lows one to perform the tomography passively, and show 
a specific method in which the entanglement of a string 
of photonic cluster state can be measured directly for up 
to 6, 20, 80 photons for the same emission/collection 
parameter values mentioned above. More importantly, 
we will show that from the same measurement data it 
is possible to indirectly infer entanglement between sev- 
eral hundreds of photons, even with standard Pauli noise 
rates per photon at the several percent level. 

The trick will be to devise a scheme that does not re- 
quire a large number of successful detections in a row. 
Fig. [1] shows our proposed experimental setup for the to- 
mography. The emitted photons pass through two beam 
splitters, randomly ending up in one of three detectors. 
One beam goes into a detector measuring circular polar- 
ization (eigenvectors denoted as |L), \R)), one in hori- 
zontal or vertical polarization (|7J),|1^)), and the third in 
diagonal polarization {\D),\D)). These three polariza- 
tion bases will be denoted as the Z basis. X basis and Y 
basis respectively. The translation from one polarization 
basis to another can be done by quarter-wave plates and 
half-wave plates, as depicted in Fig. [TJ 

The experiment now proceeds as follows: photons are 
emitted from the cluster state machine gun, which can be 
either a quantum dot in a microcavity [l| or any other 
similar setup 043 which emits photons entangled in a 
linear cluster state. These photons then go one by one 
through the beam splitters, which act classically and just 
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randomly direct each photon towards a detector. Denot- 
ing hy X,Y and Z a photon measured in the correspond- 
ing basis, and by __ a photon which was lost because 
of collection efficiency, a typical measurement result will 
look like a random string oi'X\'Y^ , 'Z' and together 
with the measured eigenvalues: 
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In some of what follows measurements in only two bases 
are needed which can be achieved with one less beam- 
splitter. 

In order to show that this passive measurement setup 
suffices to give us all the information needed for mea- 
suring the entanglement length we make the following 
assumptions: First, we presume the state of the emit- 
ted photons is translationally invariant; i.e. the state of 
photons i, . . . ,i + I is identical to the state of photons 
j, . . . , j + I, for every i, j, I. This amounts to assuming 
a continuous firing of the machine gun for long time pe- 
riods is possible, but does not amount to assuming that 
the error mechanisms are slow over such periods. The 
second assumption is that the errors of the state are in- 
dependent Pauli errors on one or perhaps two adjacent 
photons. In the proposed scheme for the pulsed optical 
cluster state as well as in most other setups, both 
these assumptions hold - they are basically a property of 
the cluster state being created. 

The localizable entanglement (LE) between two parti- 
cles is a measure of the maximal possible bipartite en- 
tanglement that can be created between them by local 
operations on the other particles with which they are cor- 
related @, 01 • In perfect linear cluster states the distance 
over which the LE is non-zero is infinite. For an imper- 
fect linear cluster state output from the machine gun, 
our goal will be to determine the distance (number of in- 
between photons) over which the LE between photons is 
non-zero. We define this quantity as £^e- Therefore find- 
ing > gives an operational indication of whether two 
particles share useful entanglement. 

Imagine now that measurement of the localizable en- 
tanglement LE{k, k + I) between photon k and photon 
k + I, which are I photons apart, is attempted. Denoting 
by a measurement on all the photons except photons 
k, k + l, (the measurement can include photons which are 
not really measured, but just traced out), it must hold 
that 



LEik,k + l)>J2PsEipik+i) , 



(2) 



with s running over the outcomes of the measurement 
A4, Ps is the probability of said outcome and E{pl. ^.^j) 
is the entanglement between photons fc, k + I, which 
can be deduced from the reduced density matrix. Thus, 
this gives us a lower bound on the localizable entan- 
glement. We will be interested in measurements that 
have the same entanglement for the 2 photons k, k + I 
for all measurement outcomes, thus reducing Eq. 1^ 
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FIG. 2. (color online) The figure depicts the measurement 
sequence Mc and the cluster state stabilizers which are re- 
sponsible for the entanglement between the two unmeasured 
photons k, k + l for 1 = 8. The Pauli matrices depicted in red 
are the stabilizers of the photons k, k + I after the measure- 
ment sequence Mc- 



to LE{k,k + I) > E{pk,k+i)- These measurement se- 
quences do not change depending on previous measure- 
ment outcomes, and can therefore be measured passively 
(see Ref. Q for a more detailed discussion). Also, be- 
cause the cluster state has translational invariance (its 
stabilzer group is translationally invariant), we can de- 
note = Pk.k+i, for all values of /, with p\2 now being 
the reduced density matrix for two photons distance / 
apart. Eq.([2]) now becomes 



LE{k,k + l) > E{p[2). 



(3) 



In order to calculate the entanglement i?(pi2), since the 
reduced density matrix will typically be mixed (as we 
have traced over some of the other photons plus there is 
possibly noise) , we will use the entanglement of formation 
M' fioj for which an explicit formula is available. 

The question we now face is: what measurement se- 
quence Ml will maximize Eq. [31 and give us the best lower 
bound on the LE? Since the experimental state will be 
a cluster state with additional errors, let us first explain 
how to get the best lower bound on the LE for the perfect 
cluster state. 

For the ideal cluster state [ll| , we use the special prop- 
erties it possesses from being a stabilizer state [l^ , with 
the stabilizers Ki = Zi-iX^Zi+i. A set of Pauli mea- 
surements denoted Mc is shown in Fig. [5] that has the 
property it collapses qubits k and k + I into a maximally 
entangled state (similar to the measurement set intro- 
duced in Ref. [13 )• To see why, consider a pattern of 
products of stabilizers of the form KiK2K4K5KiKs.... 
This yields a sequence of Pauli operators depicted as Fi 
in Fig. [5J If we shift this sequence one qubit to the left 
(recall the state is translationally invariant) we obtain 
the second form of Fi depicted. Consider also a pattern 
of products of stabilizers of the form K1K3K4KQKJ..., 
which yields a pattern F2, also shown in the figure. Both 
Fi and F2 must yield a value -1-1 when measured on the 
cluster state, since they are products of stabilizers, and 
both can be found as sub-sequences of A^c- From this 
one deduces that the resultant state on qubits k and k + 1 
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in the figure are stabilized by {YZ, ZY, XX}, which are 
tlie stabilizers of a maximally entangled state. 

Even when the relevant state is not the perfect cluster 
state we will show the measurement sequence Mc gives 
a very good lower bound on H- In this case, the 
entanglement length will be finite, but can still be large. 

There arc two parts to our proposal for measuring 
S^E- The first involves a direct determination of (a lower 
bound on) E{p{2) for any given I by finding appropriate 
sequences of measurements in the data. The second idea 
is to use a more detailed analysis of how the expecta- 
tion values of certain sequences in the data are changing 
with the length of the sequence to determine the rate 
of various errors on the cluster state. If the model fits 
(which can be verified) then very large values of can 
then be inferred. We refer to these as the direct and in- 
direct methods for inferring and begin by analysing 
the direct method. 

Key to our proposal is that instead of looking for the 
full measurement sequence A4c, as discussed above we 
can measure (Fi) and (r2). That is, taking the very long 
string of measurements in random bases outputted from 
the experimental setup of Fig. [U one can expect to some- 
times find instances of the sequence Fi or F2 from which 
these expectations can be determined. Both of these se- 
quences allow for certain qubit locations to remain un- 
determined. That is, such photons can be lost without 
detriment - we require less photons "in a row" to be col- 
lected. As is clear from Fig. [2J (Fi) gives us the values 
for Tr [pi^iZ ® Y)] = pzY and Tr [pi^iY ® Z)] = pyz, 
while (F2) - Tr [pi^iX ® X)] = pxx- 

With these values in hand we can find a lower bound 
for £'(^'12)- In fact E{p[2) > E{p\2) where we define 

p'i2 = ^{I + PYzY ® Z + pzyZ ®Y + pxxX ® X). (4) 

That is, although the actual state of the two qubits will 
(since it is imperfect cluster state) contain other terms 
in the density matrix, these can only increase the possi- 
ble entanglement of the state. To see this note that an 
LOCC protocol (which cannot increase entanglement) in 
which the three operators ZY, Y Z and XX arc each ran- 
domly applied to the two qubits always brings p\2 to the 
form of p5^2- Note also that not only are we not having 
to perform full tomography on the I qubits to demon- 
strate their entanglement, we are not even performing 
full tomography on qubits k,k + I. 

We now calculate the number of photons needed for 
a single instance of a useful measurement sequence of 
length I. Both Fi, r2 have a preferred measurement 
direction Y, in which the number of measured photons 
increases with the length of the sequence. Denote by pp 
the probability that the beam splitter passes the photon 
to the F-direction detector and by pd the probability of 
photon detection. Since the probability of getting even a 
start of a string is small, we can assume that the mean 
number of photons needed for one instance of the string is 
1 /pi , with pi being the probability of getting a sequence 
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FIG. 3. The maximum instance length for each one of the se- 
quences Fi , r2 which will be measured by the direct method 
in the experimental setup (depicted in Fig. [1]) as a function of 
the photodetection+coUection probability. The rate of emit- 
ted photons is taken as 1 photon every nanosecond, and the 
measurement time is 10 seconds. 



of length I. We get that 



"m n,, (I Pp 
Pi=Pd Pp 



(5) 



with Um being the total number of photons measured, Up 
the number of photons measured in the preferred direc- 
tion (Y) and a is the number of detectors in the experi- 
mental setup not in the preferred direction (a = 1 for Fi, 
and a = 2 for F2). Now, for each sequence, one can find 
the optimal value of Pp = by derivation. Plugging 
this back into Eq. [5] gives 




for / > 2. 

Therefore, since the number of photons needed for a 
sequence F is the corresponding 1/pf , the amount of time 
needed for one instance of length I is Te,„/pf , with T^ra 
being the the average time for the emission of one photon 
(and is proportional to T^ycie for the pulsed cluster state 
scheme). Fig. [3] shows the maximum instance length for 
each one of the sequences Fi, F2 which will be measured 
in the experimental setup (depicted in Fig. [IJ in a rea- 
sonable time (10 seconds) as a function of the photode- 
tection+coUection probability. This assumes T^^a = ^ ns 
for the quantum dot (or the emission rate for other ex- 
perimental emitters). We see that for reasonable pho- 
todetection-l-collection probabilities, sequences of length 
greater than 15 can be measured directly. 

The direct method above gives first hand information 
about the LE for the measured sequence lengths. We now 
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FIG. 4. The localizable entanglement length as a function 
of the error probability pzz, for two different values of Pauli 
error probability Pa — Q, Pcr ~ 0.2%. 

describe an indirect method for extrapolating from the 
same random sequence of measurements a much higher 
lower bound for 

To do this we need to consider the possible errors on 
the linear cluster state. The two scenarios which will be 
investigated are the case of single random Pauli errors, 
and the case of Zi ® Z^+i errors, both of which are rele- 
vant for the photonic cluster state strings In fact, as 
was shown there, single qubit Pauli errors in the emitter 
transfer into being Zi ® Zi+i errors on the output pho- 
tons. We should therefore expect these pairwise Pauli 
errors to be the most relevant. Once a photon is fired 
out then there is no mechanism to generate single Pauli 
errors on it (there is no populated bath it couples to at 
optical frequencies). However, given the variety of sys- 
tems being considered to build a machine gun it could 
well be that different schemes to quantum dots do pro- 
duce a small single Pauli error rate on the photons, and 
so we consider this case as well. 

For single Pauli errors (i.e. each error acts as applying 
a Pauli operator on a single site), with the error proba- 
bility po-, let us understand what value one should expect 
for a measurement sequence F of length I. Assuming a 
uniform probability for the three different Pauli errors, 
for each of the measured photons two out of the three 
Pauli operators will cause F to flip sign. This is true for 
any odd number of such errors. However, an even num- 
ber of such errors on the measured photons will cancel 
out. Using the binomial theorem we get that 

(r) = (i-^P.)"", (6) 

with Um being the total number of photons measured in 
the sequence F. Therefore, in this case, the less photons 
needed to be measured for a sequence of length I, the 



higher the expectation value. Thus, all the expectation 
values for the different measurement sequences will de- 
crease as I increases, but at different rates corresponding 
to the number of measured photons for each measure- 
ment sequence. Now, from the value for (F) for different 
values of / calculated from the experimental measurement 
string, it is possible to extract Pa (and see by the fit to 
the experimental data if this is the right model to use). 
Given p^r and the values of Fi, F2, it is possible to ex- 
tract the localizable entanglement for any length I, and 
find the value of ^b. 

For the case of Z ® Z errors, again assuming uncor- 
related errors with a uniform error probability pzzi we 
get an expectation value for each of the sequences that 
decreases exponentially with the length I: 

(Fi) = {l-2pzzy'\ 
(F2) = (l-2pzzY' ■ 

As for the single Pauli error case, we can use sequences 
of varying length / to infer the error rates and then lower 
bound the localizable entanglement length ^g. Fig. |4] de- 
picts the bound we can obtain as a function of the error 
probability pzz^ for two different values of single qubit 
Pauli error probability p^ = 0, p^ = 0.2%. It is im- 
portant to note that this indirect method yields a much 
tighter lower bound than that which would be achieved 
by bounding S^e directly by measuring only the cluster 
state stabilizers Ki = Zi-iXiZi+i. We see that even for 
high error rates it will be possible to confirm the gener- 
ation of long strings of entangled photons. 

We conclude with a few observations. There are other 
possible measurement sequences to consider. For in- 
stance one can use stabilizers which correspond to mea- 
suring photons k+l^k -1-2, fc -|- Z — 1 in the X basis , and 
the rest in the Z basis [IJ, [la] . These sequences actually 
give a tighter lower bound than Fi, F2 in some cases. 
However, these sequences require measurements of more 
photons, so shorter sequences will be found in the mea- 
surement data at finite collection efficiencies. We make 
no claim as to the optimality of the measurement se- 
quences we have chosen to illustrate our key ideas above. 

Finally, one can also use the experimental measure- 
ment string outputted from the experimental setup for 
other purposes: reconstructing the pure quantum state 
with the highest fidelity with the experimental state [l^ , 
or finding the fidelity of the experimental state with the 
ideal cluster state [13, [3 • Determining the efficiencies of 
these procedures as applied to the photon machine gun 
is an open problem. 
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